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Abstract 

In this paper, we study a 3D compact U(l) lattice gauge theory with a variety of nonlocal interactions 
that simulates the effects of gapless/gapful matter fields. We restrict the nonlocal interactions among 
gauge variables only to those along the temporal direction and adjust their coupling constants 
optimally to simulate the isotropic nonlocal couplings of the original model. This theory is quite 
important to investigate the phase structures of QED3 and strongly-correlated electron systems like 
the 2D quantum spin models, the fractional quantum Hall effect, the t-J model of high-temperature 
superconductivity. We perform numerical studies of this theory to find that, for a certain class 
of power-decaying couplings, there appears a second-order phase transition to the deconfinement 
phase as the gauge coupling constant is decreased. On the other hand, for the exponentially- 
decaying coupling, there arc no signals for second-order phase transition. These results indicate the 
possibility that introduction of sufficient number of massless matter fields destabilizes the permanent 
confinement in the 3D compact U(l) pure gauge theory due to instantons. 
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1 Introduction 



Gauge theories and their associated concepts have played important roles not only in elementary 
particle physics but also in condensed matter physics. For example, the conventional superconduct- 
ing phase transition is characterized as a change of the gauge dynamics of the electromagnetic U(l) 
gauge symmetry. Also, for a variety of strongly-correlated electron systems, it has been recognized 
that their phase structures and the properties of low-energy excitations are naturally described by 
using the terminology of gauge theories[l]. As such low-energy "quasi-particles" , the composite 
fermions/bosons in the fractional quantum Hall states[2] and the holons and spinous in the t-J 
model of high- Tc cuprates[3] have been proposed. It was argued[4] that their imconventional prop- 
erties like fractionality and their existence itself may be explained by a confinement-deconfinement 
phenomenon of the gauge dynamics of the effective gauge theories derived from the original models. 
This interesting idea is still controversial [5], but certainly warrants further investigation. 

Most of the problems in gauge-theoretical studies on the strongly-correlated electron systems re- 
duces to studying the phase structures of gauge theories coupled to gapless relativistic/nonrelativistic 
matter fields. In the elementary particle physics, it is generally believed that the phase structures 
of gauge systems coupled to matter fields arc not easy to study analytically because introduction 
of matter fields results in the lack of simple order parameters such as Wilson loops for pure gauge 
systems[6]. Furthermore, inclusion of gapless fermions make numerical simulations a difficult task 
since one must face quite nonlocal interactions generated by integrating over fermion variables. 
Sometimes couplings to gapless matter fields change the universality class of the gauge system under 
consideration from that of the pure gauge model. A good example is the four-dimensional (4D) 
QCD coupled with light quarks in which the number of light quarks strongly influences its phase 
structure [7]. 

In this paper, we address this problem of gauge dynamics of coupled systems. Specifically, we are 
interested in the U(l) lattice gauge theory (LGT) with gapless/gapful and relativistic/nonrelativistic 
matter fields in three dimensions (two spatial dimensions at zero temperature). This theory of course 
covers the important model QED3[8]. It appears also as a main part of the effective gauge theory 
of the strongly-correlated electron systems mentioned above, and so plays an important role in 
studying these systems[4, 9]. In the ordinary 3D compact U(l) pure gauge system (i.e., without 
matter fields) with local interactions, it is established that only the confinement phase is realized 
because of instanton condensation[ll]. For the case with additional massless (gapless) matter fields 
coupled to the U(l) gauge field, recent studies give controversial results on the possibility of a 
deconfinement phase; it is supported in Ref.[4, 12] whereas it is denied in Ref.[13]. 
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Our approach to this problem is by (i) introducing an effective theory of the original theory and 
(ii) studying its phase structure numerically. As the effective theory we use a 3D U(l) pure LGT 
with nonlocal interactions among gauge variables. These nonlocal interactions are along the temporal 
direction and mimic the effect of matter fields. We consider exponentially-decaying interactions for 
massive matter fields (i.e., fields with gaps) and power- decaying interactions for massless (gapless) 
fields. We shall see that certain cases of power-decaying interactions exhibit second-order phase 
transitions which separate the confinement phase and the deconfinement phase[10]. The existence 
of the deconfinement phase in the effective theory indicates that it is realized in the original model if 
the number of gapless matter fields is sufficiently large. This result is in agreement with the results 
of Ref.[4, 12]. 

The rest of the paper is organized as follows. In section 2, we briefly survey the effect of matter 
fields in several aspects. In section 3, the original model and its effective nonlocal gauge model 
are explained. Section 4 is devoted for numerical calculations. We calculate the internal energy, 
specific heat, expectation values of Polyakov loops, Wilson loops, and density of instantons. All 
these quantities indicate a second-order confinement-deconfinement phase transition (CDPT) for 
the gauge models with sufficiently long-range correlations among the gauge variables. In section 5, 
we study tractable low- dimensional spin models, which are obtained by simple reduction of the gauge 
degrees of freedom in the nonlocal effective model. Then we obtain an intuitive picture of the CDPT 
of the present long-range U(l) gauge theories. Section 6 is devoted for conclusion. In Appendix, 
the effective nonlocal model is studied by the low- and high-temperature expansions. The analytic 
expressions of the internal energy and the specific heat are in good agreement with the numerical 
calculations in Sect. 4. 

2 Effects of Matter Fields on Gauge Dynamics 

In this section, we briefly review the effects of matter fields upon the U(l) gauge dynamics in several 
aspects. 

2.1 Weeik-coupling regime in noncompact /compact U(l) gauge theory 

One may generally expect that inclusion of massless matter fields to a system of gauge field may 
drastically change the gauge dynamics at long wavelengths. For the case that relativistic and massless 
matter fields are coupled to a U(l) gauge field Ai^{x) {fi = 0,1,2) in 3D, the one-loop radiative 
correction of matter fields at weak gauge couplings generates the following nonlocal term in the 
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effective action of ^^(x); 

AA(xe^ [ d?x [ d^yJ2Fi.Ax)^^^F^,{y), F^, = d^A, - d^A^. (2.1) 

It is obvious that the above term strongly suppresses fluctuations of ^^^(a;) at long distances due 
to the factor \x — The above nonlocal terms are leading at low energies and momenta «/its 

eff'ective coupling constant does not vanish at the infrared limit. Then the potential energy V{r) 
between two charges separated by distance r changes drastically, 

y(r) oclogr y(r) oc -. (2.2) 
r 

In the compact U(l) gauge theory, topologically nontrivial excitations (e.g. instantons) of j4^(a;) 
appear whose effect at small gauge coupling e can be estimated by replacing the field strength 
Ffivix) in Eq.(2.1) by Ffj_i,{x) — 27rn^,^(a;), where n^,^ is an integer field whose rotation measures 
the instanton number (density) p{x) = e^v\dij.nv\{x). In the gauge model with the usual Maxwell 
term, the potential energy between a pair of instantons at distance r is Coulombic, Vins(?') ~ 
However, the long-range action (2.1) modifies Vins(?') to a long-range one[12], 

T4ns(r-) (X-^ 14ns (0 oc log T. (2.3) 

r 

Recently, it was argued that a gas of charged particles with the long-range interaction like logr 
exhibits a phase transition between a dilute gas of dipoles and a plasma[14]. This result implies 
that instantons in the long-range gauge theories form dipole pairs and do not condense in the weak- 
coupling region. Since the confinement phase of gauge dynamics requires a condensation of isolated 
instantons (a plasma phase), this result leads to a deconfinement phase of gauge dynamics there. 



2.2 The 3D CP^"^ model on the criticality 

As an example of the nonperturbative effect of matter fields, let us consider the CP^~^ model in a 
3D continuum. The action of the model is given as 



j Sx[-^\{d^ + iA^)z\^ + a[\z\''-l)\, (2.4) 



where Za{x) (a = 1, 2, • • • , TV) is the CP^~-^ field satisfying the local constraint, Y^^=i \^a{x)\'^ = 1 
for each x via the Lagrange multiplier field a{x), and A^{x) is the auxiliary U(l) gauge field. 

At large N, the expansion is reliable and predicts a second-order phase transition at the 
critical coupling f = fc [15]. For / > fc, the model is in the disordered-confinement phase in which 
{<y{x)) ^ 0, whereas for / < /c, the model is in the ordered-Higgs phase in which {<j{x)) = and the 



4 



QpAT-i ggjjj Yias a nonvanishing expectation value like {zn{x)) = vq ^ 0- As a result, the gauge field 
Af^ acquires a finite mass (cx vq) by the Anderson-Higgs mechanism, and the low-energy excitations 
are gapless Za{x) {a ^ N) fields. 

Recently, considerable interests have been paid on the question how the gauge field behaves just 
at the critical point f = /c[16]. In the leading order of the 1/N expansion, it is shown that all 
the components Za{x) (a = 1, 2, • • • , A'') are massless and the gauge field acquires the nonlocal 
"kinetic term" of Eq.(2.1). This implies that, at the critical point, the nonperturbative fluctuations 
of like instantons are suppressed, so the gauge dynamics at / = /c is in the Coulomb phase with 
the potential V{r) ~ 1/r as Eq.(2.2). 

Then we have numerically studied the 3D CP^+ U(l) LGT from the above point of view[17]. 
The numerical results of the CP^ model show a second-order transition and the conflnement phase 
is realized at the critical point. However, calculations of CP^~^ model with A' = 3,4, 5 show that 
the topologically nontrivial configurations are suppressed more as A' increases. We expect that there 
exists a critical value of A' = A^c and the deconfinement phase is realized at the critical point for 
Nc < N in compatible with the large-A' analysis. We note that a similar phenomenon of inducing 
a deconfinement phase by a plenty of massless matter fields has been established in lattice QCD in 
4D. For a suflaciently large number of light flavors Nf > 7, the model stays in the deconfinement 
phase even the pure gauge term is missing [7]. 

2.3 Nonrelativistic fermions in strongly-correlated electron systems 

Nonrelativistic fermions are distinguished from relativistic fermions by the properties; (i) they prop- 
agate only in the positive direction of the imaginary time in path-integral formulation, and (ii) 
they form a Fermi surface (line). In strongly-correlated electron systems, one faces nonrelativistic 
fermions not only in the original models but also in their effective gauge models. In studying the frac- 
tionalization phenomena of electrons like the charge-spin separation (CSS) [3] in high-temperature 
superconductivity and the particle-flux separation (PFS) in fractional quantum Hall systems [2] , we 
regard an electron Cx (we suppress the spin index for simplicity) at the site a; of a 2D spatial lattice 
as a composite of a fermion A^ and a boson [4] as 



To assure the correct physical space composed of Cx, one imposes the local constraint for the physical 
states |phys) as 



Cx = AiB, 



(2.5) 



[AIAx + BIBx - 1) Iphys) = for each x. 



(2.6) 
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The hopping term of electrons (e.g., the t-tevm of the t-J model) may be rewritten via a "decoupling" 



as 



iBl^,W,iB,+Al^,wl,A, + ii.c.)-\W,i\^ + ---. (2.7) 

Wxi is an aiixiliary complex field defined on the link (a;, x + i) where i = 1, 2 is the direction index 
of the lattice. Its phase degree of freedom Uxi = Wxi/\Wxi\ behaves as a spatial component of a 
compact U(l) gauge field Uxi[= exp(i^xi)]- Uxi represents the binding force of the constituents Ax 
and Bx- In the path-integral formalism, the partition function Trc = exp[— /3if(C)] (/? = T~^) has 
the following representation: 

Z = f HdAxiTjdAxiTjdBxir) H dC/^^(T)exp(A), 

•' x,T /i=0,l,2 

A = dry^l - iOxo - Ax{do - iOxo + IJ'a)Ax - Bx{do - iOxo + Mb)-Bx 

•^0 X ^ 

+t {A^+iUxi{T)Bx + H.C.) + Ant] , (2.8) 

i 

where r (g [0,/?]) is the contimium imaginary time, Ax{t) is a Grassmann number expressing 
fcrmions, ^.a(b) is the chemical potential, and t is the hopping amplitude. The field Oxo is a Lagrange 
multiplier field to enforce the constraint (2.6). It may be viewed as the time-component of gauge 
field. In fact, in the discrete-time formulation on the 3D lattice it appears as the exponent of gauge 
variable in the r-direction, Uxo = c'X-'p{i9 xo) ■ 

When the gauge dynamics of Uxi is realized in the confinement phase, the constituents Ax and Bx 
are bound within electrons and the relevant low-energy quasi-particles are described by the electron 
operators Cx themselves. On the other hand, if the gauge dynamics is realized in the deconfinement 
phase, the binding force is weak and these constituents are no more bound in each electron but 
dissociate each other. This is just the the fractionalization(separation) phenomena of electrons. 

In Ref. [4] we have studied the system of Eq.(2.8) at finite temperatiues (T) by the hopping 
expansion in the temporal gauge 6xa = 0. After integrating over Ax and Bx in powers of i, one 
obtains the effective interactions. Up to 0{t^), by restoring the temporal components Uxq-, one gets 
the following effective interaction among Ux^, 

AA oc 5(1 - (5) • V / dT dr' Vx^ i{T, r') + c.c, 

x,i 

^xi,i(T,r') = Ux^,r',i exp{i J c?t"[6»^^_^j o(r") - 6»a;^,o(T")]) C/xi,r,i, (2.9) 
where 5 = (A^Ax) is the concentration of fermions. The corresponding terms are illustrated in Fig.l. 
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Figure 1: Ilhistration of the nonlocal interactions in electron- fractionalization plienoniena. Px^@ = 
exp(i J drOxo) is the Polyakov line at the transverse coordinate x± = (.ti,.X2), Qx±.i{'ri,T2) is given 
by the product of propagators {Bx^{ti)BI^{t2)) and {Ax^{t2)AI.^{ti)) . The product of these Px^^ 
and Qx±,i{Ti, T2) results an effective interaction V^^,i(Ti, T2) along a closed loop, which has the same 
form as a vacuum polarization of relativistic matter fields. 



The interactions among Uxfiir) in A A are quite nonlocal in the r-direction, and favors the 
ordered configurations of Uxfi, so the deconfinement phase. In the previous papers[4], we argued 
that this is the essence of the mechanism of fractionalization phenomena like CSS and PFS. By 
mapping these gauge models approximately to a spin model, we concluded that the U(l) gauge 
dynamics is realized in the deconfinement phase at the low-T region below certain critical line Tc{S). 
Thus the possible deconfienment phase in the model similar to Eq.(2.9) is to support the electron 
fractionalization phenomena. 



2.4 Chern-Simons term by Dirac fermions 

On considering the effects of matter fields upon 3D gauge dynamics, there is an important difference 
between fermionic and bosonic matte fields. That is, relativistic fermions have a possibility to 
generate the Chern-Simons (CS) term A^cs via radiative corrections: 

AAcs = cjd^x ei^^pAi^d^Ap, (2.10) 

which violates the parity symmetry. In particular, as the mass term mo'tptl' of 3D two-component 
spinor Dirac fermion field ip violates the parity invariance, it generates the CS term with the coef- 
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ficient c oc sgn(m£)) [18]. On the other hand, scalar fields do not renormalize the coefficient of the 
Chern-Simons term. In the perturbation theory, the CS term is the leading term at long distances 
and low energies. So it may change the phase structure of the gauge system, in particular that of 
the compact gauge systems. One can intuitively expect appearance of a deconfinement phase since 
the CS term suppreses fluctuations of the gauge field. 

However in most of the effective gauge theories of strongly-correlated electron systems, the CS 
term is not spontaneously generated, because the most of these systems including the t-J model 
preserves the parity invariance. For example, in the flux state of the Heisenberg antiferromagnetic 
spin model and the t-J model in the slave fermion representation, relativistic Dirac fermions appear 
as low-energy excitations, but the parity invariance is preserved in the effective theory because the 
lattice fermions appear in doublets with opposite signatures of masses [9], which correspond to four- 
component spinor Dirac field in the continuum. Then the CS coefficient cancels with each other in 
the radiative correction as c cx X^^^^j.^ sgn(m£)) = 0. 

Though nonperturbative investigation of the Chern-Simons gauge theory is very important by 
itself^ , we shall focus on parity-invariant lattice gauge theory in the rest of discussions in this paper. 

3 Nonlocal U(l) Lattice Gauge Theory 

In the previous section, we have seen various approaches to study the effect of matter fields upon 
U(l) gauge dynamics. To confirm these results, one must examine the validity of the approximations 
employed there. In particular, to check the validity of the hopping expansion at T = is quite 
important for studies on strongly-correlated electron systems. 

Keeping this problem in mind, we start this section with a U(l) LGT coupled with matter fields 
and introduce its effective nonlocal LGT. Let us consider a 3D cubic lattice (i.e., a 2D lattice with a 
discrete imaginary time). The gauge field Ux^ {n = 0, 1, 2) is defined on the link {x, x + fi) between 
the pair of nearest-neighbor sites x and x + fi. The partition function Z is given by the following 
functional integral, 

Z = Y]_d(j)xd(l)a:YldUxij.e:icp{A), 

X Xfl 
A = -^4>xTxy{U)(t)y+Au, 

x,y 

= q ^ iUxuUx+i>,nUx+fi,Mxn+ C.C.), (3.1) 

^For example, phase structure of SU(N) Maxwell-CS gauge theory has been studied via frustrated Heisenberg spin 
model without parity invariance, which is a low-energy effective model of strongly-coupled SU(N) gauge theory of 
fermions [19]. Existence of a deconfinement phase transition is suggested there. 
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where x = {xo,xi,X2) is the site-index of the 3D lattice of the size V = N0N1N2 with the periodic 
boundary condition, /x is the imaginary-time index {fj, = 0) and spatial direction indices (/U = 1,2), 
is the matter field on x, Ux^ = exp{i9xij,) (— tt < Ox^ < tt) is the U{1) gauge variable on the link 
{x, X + jl), and q is inverse gauge coupling constant. Txy{U) represents the local minimal couplings 
of (j)x to Uxn- For example, for a bosonic matter field, 

Y,4>xTxy{U)<t>y = tY,[4>x+f,Ux^<t>x+'^-c]+ M^Y,^^'t>x, M^ = Q + m^, (3.2) 

x,y x^n X 

where w? is the mass in unit of the lattice spacing and 6 in is the number of links emanating 
from each site. 

After integrating over the matter field (j)x-i effective gauge model is obtained, which includes all 
contributions from (j)x to the gauge dynamics. 



/n 



dUx^, exp \J Tr log Txy{U) + Au\ , (3.3) 

where / is a parameter counting the statistics and internal degrees of freedom of (j)x- Due to the 
(Tr \ogV xyiU)) term, the effective gauge theory becomes nolocal. For relativistic matter fields, a 
formal expression of the eflFective gauge theory action is obtained by the hopping expansion and it 
is expanded as a sum over all the closed random walks TZ (loops including backtrackings) on the 3D 
lattice, which represent world lines of particles and antiparticles as 

TT\ogTxy{U)=^jj— n Ux^- (3-4) 

n ^ ' {xn)eTZ 

L[JZ] is the length of TZ, and 7 = (6 + m^)~^ is the hopping parameter. There are many difFercnt 
random walks that have the same shape of a closed loop on the lattice. Each random walk in such 
a family may have a different starting point and/or backtrackings. This degeneracy cancels out the 
denominator L[TZ\ in Eq.(3.4). For the constant gauge-field configuration Ux^ = 1, the expansion in 
(3.4) is logarithmically divergent ^ logm as m ^ due to the lowest-energy zero-momentum mode. 

Below we shall study a slightly more tractable model than that given by Eq.(3.3). It is suggested 
by the hopping expansion (3.4), and obtained by retaining only the rectangular loops extending in 
the r-direction in the loop sum and choosing their expansion coefficients as follows; 



Zt = Y[dUx^exp{AT), 

2 .Yo 

At = Cr{Vx,i,T + Vx,i,r) + As, 

x i=l T=l 
r-1 

Vx,i,T = Ux+tQ,i ri [^x-|-fe6,0^x-|- 



|-i-|-feO,0 

k=0 " 



Uxit 



9 



As = 5AXl(f^^2?7^+2,iJ7^+i,2C^xi + c.c.), 



(3.5) 



where g is the (inverse) gauge couphng eonstant, and V^a.t is the product of U^fj^ along the rectangular 
{x,x + i,x + i -\- rO, x + rO) of size (1 x r) in the (i — 0) plane. See Fig. 2. 



=0 



C2 



X x-\-i 



Figure 2: Illustration of the nonlocal interaction in the action (3.5). Each rectangle with thick lines 
represents Vx,i,T- 



In As, we have retained only the single-plaquctte coupling with the coefficient {gX). For the 
nonlocal coupling constant c^, we consider the following three cases; 



T~", power — law decay (PD— a) [a = 1, 2, 3), 
e"™"^, exponential decay (ED), 
1, no decay (ND). 



(3.6) 



The power a = 1 in the PD model in (3.6) reflects the effect of the relativistic massless excitations 
without dimensional parameters. In fact, this Ct generates a logarithmically divergent action for 
Uxfi = 1 explained below Eq.(3.4) as one can sec from the relation, cxp(— mT)r~"'^ ~ log(l/m). 
The action for m = is then proportional to t^^ ~ logiVo for finite A^o- On the other hand, the 
ED model contains the parameter m with mass dimension and simulates the case of massive matter 
fields. (We used Eq.(3.6) instead of exp(— mT)/r to make the comparison with the PD case more 
definitive.) The ND model corresponds to gauge model coupled to nonrelativistic fermions with a 
Fermi surface (or Fermi line) as it is seen from Eq.(2.9). 
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4 Numerical Results 



In this section, wo report our numerical calculations of the internal energy, the specific heat, the 
Polyakov lines, etc., and determine the phase structure of the models. Most of our interest concerns 
a possible CDPT in the 3D compact U(l) gauge models with the nonlocal interactions. 

In the MC simulations, we consider the isotropic lattice, = N (p, = 0, 1, 2), with the periodic 
boundary condition up to iV = 32, where the limit N ^ oo corresponds to the system on a 2D 
spatial lattice at T = 0. For the mass of the ED model, we set m = 1. For the spatial coupling A 
scaled by g, we consider the two typical cases A = (i.e., no spatial coupling) and A = 1. 

4.1 Internal energy and specific heat 

First we calculate the following "internal energy" E and the "specific heat" C per site; 

Zr = [ [dU] exp(Ar) = exp{-FV), 



Here we note that the definition of C (4.1) is the response of E under the variation of g and is different 
from the conventional specific heat that measures the response under the variation of temperature 
itself. The latter contains extra terms associated with the change of Nq. Because these terms behave 
less singular than the variance of E, they are irrelevant in searching for phase transitions. 

In Fig. 3, we present E and C for A = 1 vs the nonlocal coupling g. In Appendix we calculate E 
and C by the high-temperature expansion (HTE) for small g and by the low-temperature expansion 
(LTE) for large g. The MC results of E and C are consistent with these expansions. For comparison 
with HTE (e.g., for < 51 < 0.1 in the PD-1 model {a = 1)), see Fig.15 in Appendix. For 
comparison with LTE, the leading result of LTE is shown by the straight lines in Fig.3(a,b) and 



First, let us see the PD-1 model (a = 1) in detail. E of Fig. 3(a) connects the results of HTE 
and LTE. Since the LTE is an expansion around Vx,i,T = 1, this behavior implies y^.i.r ~ 1 for large 
g. C of Fig. 3(c) shows that its peak develops as the system size A'' increases. These two points 
indicate that the PD-1 model exhibits a second-order phase transition separating the disordered 
(confinement) phase and the ordered (deconfinement) phase at <? = (?c — 0.17 which is determined 
from the data of A'' = 24. The existence and the nature of the phase transition will be confirmed 
later by the measurement of the Polyakov lines and the instanton density as we show in the following 
subsections. 



E 



C 




(4.1) 



C = 1 in Fig.3(c,d). 
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Figure 3: Internal energy E and its fluctuation C of the action with A = 1 vs non-local coupling g for 
N = 8(4), 16(*), 24(B); (a,c) PD-1 model, (b,d) ED model. The solid lines in (a) and (b) are the 
leading order in the large-g expansion (LTE). In the PD-1 model, strong N dependence is observed 
in E at large g and in the developing peak of C. They indicate a second-order phase transition in the 
PD-1 model. In the ED model, on the other hand, C does not develop and therefore the observed 
peak indicates not a phase transition but a crossover. 
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We notice that the location of the peak in C of Fig. 3(c) shifts to smaller g as the system size N 
increases in the direction opposite to the usual second-order phase transitions. This behavior reflects 
the fact that the couplings among increases effectively as A'' increases due to the additional 
terms in the summation over A^o in the action even if one fixes the overall constant g. This is the 
characteristic nature of nonlocal interactions in strong contrast to local interactions. 

On the contrary, in the ED model of Fig. 3(d), the peak of C does not develop as A'' increases, 
showing no signals of a second-order transition. It may have a higher-order transition or just a 
crossover. Similar behavior of C is observed in the ordinary U(l) gauge systems with local actions 
which have only the confinement phase. The physical meaning of the above "crossover" in the ED 
model shall become clear by studying instantons in Sect. 4. 3. 

It is quite interesting to see whether the data of C for N = 8, 16 and 24 in Fig. 3(c) exhibit the 
finite-size scaling behavior [20]. To this end, let us introduce a parameter 

e = (S - fl'oo)/5oo (4.2) 

where ^foo is the critical gauge coupling of the infinite system at A'' — > oo. Then let us assume that 
the correlation length ^ scales as ^ oc e~'' with a critical exponent u. We also expect that the peak 
of C diverg es as Cpeak oc e ^ as A/^ — > oo with another critical exponent cr. The finite-size scaling 
hypothesis predicts that the specific heat C{€,N) for suflaciently large A'' scales as 

C7(e,Ar) = AT'^/XiVi/'e), (4.3) 

where ^{x) is a certain scaling function. In Fig. 4, we present ^(a;) determined by using the data in 
Fig. 3(c) with u = 1.2, ajv = 0.25 and g^a = 0.11. This result indicates that the finite-size scaling 
law holds quite well. Cbonsidering the errors in the data we estimate the values of scaling parameters 
and goo as 

= 1.2 ~ 1.3, a/i/ = 0.25 ~ 0.26, g^c = 0.10 ~ 0.12. (4.4) 

The simulations of the PD-1 and ED models with A = give similar behaviors of E and C as 
the A = 1 case, preserving the above phase structure for A = 1. That is, the PD-1 model (A = 0) 
has a CDPT, while the ED model (A = 0) has no transition. 

Let us next consider the ND model. In Fig. 5, we present C for A = and 1, which show 
strong signals of a second-order phase transition. However, the value of the critical coupling for 
A = 1, 5c ~ 0.1(A/' = 8),0M5{N = 16),0m{N = 24), decreases very rapidly as the system 
size increases. This behavior is explained by the increase of efii'ective coupling explained above. 
Then one may think that 5c — * as A'' — > oo, i.e., only the deconfinement phase survives in the 
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Figure 4: Scaling function (/){x) of (4.3) obtained from the data C of the PD-1 model in Fig. 3(c) for 
= 8, 16 and 24. The result obviously shows that the finite-size scaling law holds and the observed 
phase transition in the PD-1 model is of second-order. 



ND model. This expectation is consistent with the fact that the coefficient Q2 = J2t '^^ HTE, 
C ^ 2(2(52 + A2).g^ -I- 0(.g^) in Eq.(A.8) diverges as Q2 oc iV ^ cx) for the ND model, i.e., the radius 
of convergence in the HTE is zero. On the contrary, Q2 is finite for the PD-1 model, assuring us 
gc =/= 0. This is supported also by the scaling analysis given above. 



— 






(a) 










n 
















1k n 





















(b) 


ir- 


□ □ 






: i 


r — ° 


[ A A-A-z 












tV □ 


A 

A 















0.04 0.08 0.12 g 0.04 0.08 0.12 g 



Figure 5: C vs g in the ND model with (a) A = and (b) A 1 for N 8(a), 16(H), 24(*). The 

peaks develops as the system-size increases. 



From these results for the various cases of long-range interaction, it seems that there exists a 
critical power a — ac below which the CDPT takes place. In Fig. 6 we show C of the PD model with 
a = 2 and a — 3. 

We obtain a very interesting result, i.e., for the PD-2 system (a = 2) with the nonvanishing 
spatial coupling A = 1, there exists a second-order CDPT as in the PD-1 case, whereas in the PD-3 
case {a — 3) the peak in C does not develop as N increases, hence no signals of CDPT. Furthermore, 
careful study of the PD-2 case shows that the second-order CDPT disappears for vanishing spatial 
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Figure 6: C vs g in the PD-2 and PD-3 models for N = 8(a), 16(H), 24(*). (a) PD-2(A = 1), (b) 
PD-2(A = 0), (c) PD-3(A = 1), (d) PD-3(A = 0). The peak in C develops in (a) as N increases, 
whereas it does not in (b), (c) and (d). This indicates that the CDPT starts to appear between the 
two cases, (a) PD-2(A 1) and (b) PD-2(A = 0). 



coupling A = 0. Thus we conclude that the critical power is ac = 2 and the spatial coupling A controls 
the existence of the CDPT. The above results will be confirmed by the study of the Polyakov line 
in the following subsection. 

4.2 Polyakov lines 

We have argued the possible CDPT by measuring the thermodynamic quantities like E and C. In 
order to study the CDPT in more details and further the nature of gauge dynamics in each phase, 
it is useful to calculate order parameters in gauge theory like Polyakov lines and Wilson loops. 

First, let us introduce the Polyakov lines Px^ for each spatial site x± = {xi, X2) and study their 
spatial correlations fp{x±), 

No 

Px^ = Y\. U^±,^o,0, X± = {xi,X2), 

fpixi_) = (Px^Po). (4.5) 

Since the present model (3.5) contains no long-range interactions in the spatial directions, fp{x±^) is 
expected to supply a good order parameter to detect a possible CDPT. In the deconfinement phase, 
the fluctuations of Uxo are small, which implies an order in fp{x±), i.e., we expect fp{x±) 7^ as 
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x± — > large in the deconfinement phase. 

In Fig. 7, we present fp{x±) for the PD-1 and ED models. The PD-1 model of Fig. 7(a) clearly 
exhibits an ofF-diagonal long-range order, i.e., lima;^^oo /p(a;^) 7^ for g > 0.20, whereas the ED 
model of Fig. 7(b) does not for all g's. To sec this explicitly, we plot in Figs. 7(c) and (d) the order 



parameter p = {fp{x^-^^)y^^ for the PD-1 model, where 



.MAX 



is the distance at which 



fp becomes minimum due to the periodic boundary condition, p of the PD-1 model(A = 0) starts to 
develop continuously from zero at g = gc — 0.15. The size dependence of p shows a typical behavior 
of a second-order phase transition. Thus the gauge dynamics of the PD-1 model is realized in the 
deconfinement phase for g > gc, whereas it is in the confinement phase for g < g^. In contrast, the 
ED model stays always in the confinement phase. These results including the value of are in good 
agreement with those derived from the data of E and C given in Fig. 3. 
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Figure 7: Correlations of Polyakov lines, fp{x±), vs (a) PD-1 (A — 0, N=:16) with g ~ 

0.4,0.3,0.2,0.1,0.02 from above, (b) ED (A = 0, N==16) with g = 2.5,2.0,1.5,1.0,0.5 from above. 
In (c) and (d) the order parameter p = {fp{x^^-^)y^^ vs g is plotted for the PD-1 model; (c) PD-1 
A = and (d) PD-1 A = 1. They exhibit long-range orders for g > g^ — 0.15 in the PD-1 model for 
both A. 



Let us turn to the PD-2 model with and without the spatial coupling. In Fig. 8, we show the 
result of p for A = 1 and A = 0. We observe that the model with A = 1 shows a typical behavior of 
the second-order phase transition as the system size is increased, whereas the case of A = does not. 
From this result and the observation of C in the previous subsection, we conclude that the CDPT 
exists in the PD-2 model (A = 1) whereas it disappears in the PD-2 model (A = 0). 

In Fig. 9 we present fp{x±) of the PD-3 model. It is obvious that there is no long-range order in 
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Figure 8: Off-diagonal long-range order of Polyakov lines, p vs \xj_\ for the PD-2 model with (a) 
A = 1 and (b) A = 0. In the A = 1 case, p exhibits a steeper jump for larger N, whereas it does not 
in the A = case. 

the PD-3 model both for A = 0, 1 and only the confinement phase is realized for all g. 
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Figure 9: /p(a;±) vs g for the PD-3 model with N = 16; (a) A = 1 for g = OA ^ 0.7, (b) A = for 
g = 0.56 ~ 0.78. The upper curves have larger g's. It is obvious that there is no long-range order in 
the PD-3 model regardless of the values of g and A. 



4.3 Wilson loops 

Let us turn to study of the Wilson loops. For ordinary pure and local gauge systems, the Wilson loop 
W[C] along a closed loop C on the lattice is a good order parameter to study the gauge dynamics; 
iy[C] obeys the area law in the confinement phase and the perimeter law in the deconfinement 
phase; 

T-r cxp(—aS\C]), area law, 

W[C]^{llU,,)^\ (4.6) 
C I exp(— a'L[C]), perimeter law. 
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where S[C] is the minimum area of a surface, the boundary of which is C, and a and a' are constants. 
For a (local) gauge theory containing matter fields of the fundamental charge, W[C] cannot be an 
order parameter because the matter fields generate the terms Yic Ux^i, with coefficients ~ exp(— 6L[C]) 
in the effective action. However, in the present model (3.5), the nonlocal terms are restricted only 
along the temporal direction, so it is interesting to measure W[C\ for the loops lying in the spatial 
(1-2) plane. If W[C] obeys the perimeter law, fluctuations of the spatial component of gauge field is 
small. 




Figure 10: Wilson loops [N = 32) in the 1-2 plane at large g vs L[C] or S[C]. (a) PD(A = 1,5 = 0.25), 
(b) ED(A = 1,5 = 1.5). The PD model seems to prefer the perimeter law, whereas the ED model 
prefers the area law. 



In Fig. 10, we plot VF[C]. For the PD-1 model in Fig. 10(a), the data at 5 = 0.25 seem to prefer 
the perimeter law. For the ED model in Fig. 10(b), the area law fits W[C] better than the perimeter 
law at g = 1.5; a considerably larger value than g ~ 1.0 at the peak of C. This suggests that the 
area law holds in the ED model at all g. These observations are consistent with the previous results 
on the (non)existence of the CDPT. We conclude that the Wilson loops in the spatial plane can be 
used as an order parameter of gauge dynamics in the present model. 

The result that the PD-1 case of the simplified model (3.5) exhibits the CDPT strongly suggests 
that the original model (3.3) of massless matter fields also has the deconfinement phase, because the 
isotropic distribution of the nonlocal gauge couplings in the original model should give the similar 
effect of suppression of fluctuations of Ux^i as those in the temporal direction in Eq.(3.5). This 
expectation is supported by the measurement of the spatial Wilson loop given above. 
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4.4 Instantons 

It is well known in a 3D continuum space-time that instanton (monopole) configurations of U(l) 
gauge field carry nontrivial topological numbers, and the instanton density serves as an index to 
express the disorderness of gauge field [11]. To see the details of th gauge dynamics of the present 
model and to support the conclusions obtained in the previous subsections, let us study instantons 
on the lattice. We employ the definition of the instanton density px at the site x in U(l) lattice 
gauge theories by DeGrand and Toussaint [21]. We introduce the "vector potential" 6xn as the 
exponent of Ux^ = exp{i6xti) [Oxn € (— 7r,7r)]. Then, the magnetic flux @x,ij.i^ penetrating plaquette 
{x,x + n,x + + VjX + is expressed as 

Qx,iJ.u = 0xn + 0x+n,v-0x+i,,n-0xu, (-477 < ©x,^L/ < 47r). (4.7) 

We decompose Qx,nv into its integer part 2-Knx,^v {nx,nv is an integer) and the remaining part 
^x,tj.u = Qx.ixv (mod 2-it) uniquely, 

Qx,^lv = 2nnx,^,„ + Qx,^lv, (-tt < Ox./ii. < tt). (4.8) 

Physically speaking, n^./ai/ describes the Dirac string whereas Qx,nv describes the fluctuations around 
it. The quantized instanton charge px at the cube around the site 5 = a; + (1 + 2 + 3)/2 of the dual 
lattice is defined as 

Px ~ 2 ^ ^ ^'flL'piP'X+flfl'p "^XyVp) 

~ 4^ ^ ^ £/ii/p(0a:+M,!^p ~ ^x,vp)i (4-9) 

where e^^p is the complete antisymmetric tensor. Px measures the total fiux emanating from the 
monopole(instanton) sitting at x. Roughly speaking, Px measures the strength of nonperturbative 
gauge fluctuations around x. For the local 3D U(l) compact lattice gauge theory without matter 
fields, the average density 

p^^Ed/^-D' (4-10) 

X 

is known to behave as p a exp(— eg) (c is a constant) if the instanton action eg is large; The instanton 
gas stabilizes the confinement phase for all the gauge coupling [11]. 

In Fig. 11 we present p as a function of g in the PD-1 and ED models. It decreases as g increases 
more rapidly in the PD-1 model than in the ED model. This difference in behavior is consistent with 
the result that the PD-1 model exhibits a second-order transition, while the ED model does not. 
The A coupling enhances the rate of decrease in p as one expects since the spatial coupling enhances 
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the ordered deconfinement phase. In the ED model with A = 1, p is fitted by exp(— eg) in the dilute 
(large g) region, and the smooth increase for smaller g indicates a crossover from the dilute gas of 
instantons to the dense gas, just the behavior similar to the case of pure and local lattice gauge 
theory [11, 22]. 

In Ref. [10] we presented the figure (Fig. 3) which is very similar to Fig. 11. However, the former 
was plotted by using the definition p = V^^ Sa;((l ~ '^o.pi,))i the average of occupation number of 
instantons (1 for Px ^ and for p^ = 0) per site instead of the instanton density of Eq.(4.10) itself. 
The curves in two figures are almost indistinguishable, because the configirations with > 2 are 
rather rare. 
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Figure 11: Average instanton density p vs g. (a)PD-l(A = 0, iV = 8), (b)PD-l(A = 0, 16), 
(c)PD-l(A = I, N ^ 8), (d)PD-l(A = 1, TV = 16), (e)ED(A = 0, iV = 8,16), (f)ED(A = 1, 
N ^ 8, 16). The solid curve is oc exp(— eg) and fits (f) at large g. 

In Fig. 12 we present snapshots of p^ for the PD-1 model with A = 1. Fig. 12(a) is a dense gas 
whereas Fig. 12(b) is a dilute gas. They are separated at gc — 0.20, the location of the peak of 
C for = 16. In Fig. 12(b), instantons mostly appear in dipole pairs at nearest-neighbor sites, 
Px ~ ±l,px±^ = tIj while in Fig. 12(a), they appear densely and it is hard to determine their 
partners. In both cases, the distributions px have no apparent anisotropics like column structures. 
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However, the orientations of dipoles in Fig. 12(b) are mostly 92%) in the temporal direction as 
expected from the nonlocal interactions of Eq.(3.5). 

Next, let us examine the difference of gauge dynamics in various cases of nonlocal interactions in 
detail. To this end, it is useful to measure new observables made out of gauge-field configurations; 
nonlocal instantons elongated in the temporal direction with the length t = 2,3, - ■ ■ ,N — 1. The 
density of these instantons is defined by 

t-i 

Px,t = y^Px+rO' 

T=0 

Pt = {p.,t). (4.11) 

In Fig. 13 we plot pt vs g for various cases of the exponent a and A. 

We notice that pt's in Fig.l3(a,b,c) behave quite differently from pt's in Fig.l3(d,e). For the 
former cases, all the curves of pt with different t's decrease similarly as g increases [except for t = 1 
in (c)], whereas for the latter cases, each pt decreases in a different manner. The present nonlocal 
instanton density pt is capable to distinguish the cases (a,b,c) exhibiting a CDPT and the cases 
(d,e) without transition (cross over) through its t dependence. This is consistent with our common 
understanding that the second-order phase transition is a collective phenomenon at which wild 
fluctuations of all the variables in the system in the disordered phase start to be reduced coherently 
and also scale-invariantly near the transition point. 
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Figure 13: Density of nonlocal instantons pt vs g for = 12. (a) PD-1(A = 1), (b) PD-1(A 
(c) PD-2(A = 1), (d) PD-2(A — 0), (e) PD-3(A = 1). pt at a fixed g increases as t increases. 
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5 Effective Models at Large g and A = 



In this section, wc study a ID XY spin model for the spatial gauge variables Uxi and also a 2D XY 
spin model for the tem,pora,l ones Uxq, both of which are regarded as effective models for the present 
nonlocal model at large g and A = obtained by a simple reduction of the dynamical degrees of 
freedom. Two models are complementary each other and study of them helps us to understand the 
properties of the present nonlocal gauge model studied by the numerical simulations in section 4. In 
particular, these models are capable to explain the existence of the deconfinement (ordered) phase 
in the PD-1 model at large g and the nonexistence of it in the other PD and ED models with A = 0. 

5.1 The ID XY model with spatial variables U^i 

Let us imagine the situation in which fluctuations of the temporal gauge fields Uxo are small so that 
one may replace Uxo = exp{i6xo) by its average as follows; 



This is expected for sufficiently large g. Furthermore we put A = 0, i.e., there are no direct interac- 
tions among spatial gauge variables Uxi — exp{i9xi)- Then the system is decoupled to ID subsystems 
defined at each spatial link {x,x + i). The subsystem at {x,x + i) is described by the U(l) angles 
0j{= Oxi) (we write = J = 1, • • • , A^o and suppress the suffix Xi,X2, %)■ Its energy takes the form 
of a ID XY spin model with nonlocal interactions with couplings Cr- The partition functions of the 
subsystem and the total system are given as follows; 



To study the correlation function (exp[^(^^r — ^o)]): we make the harmonic approximation for large 
g; (i) expand the cosine term up to the quadratic term, (ii) extend the range of 9j to (—00,00) and 
(iii) neglect topologically nontrivial configurations. We expect that the third approximation in the 
above is justified by the long-range ferromagnetic interactions in Eq.(5.2). Then we have 



Uxo ^ u {Q < u < 1). 



(5.1) 




Zr(ff:large,A = 0)~(ZiDXY)' 



(5.2) 




r 
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/(r) = {e^p[i{Or-eo)])=eM-^{iOr-0of)] 
= exp[-J2G-\k){l-cos{kr))], 



(5.3) 



where 9k is the Fourier-transformed variable of 9j. For the standard local coupling Cr = Sti, the 
fc-sum in the exponent of /(r) gives / dk cos{kr) / k"^ r, which implies /(r) ~ ex.p[—{gu^)~^r] 
as r ^ oo due to the severe infrared fluctuations. For the nonlocal cases we have 

1 — cos{kr) 



f{r) ~ exp 



M 



E 



exp(— Mr), 



2r 



CtT^ 



0, 

^ 

gu 



(5.4) 



glii(l-u^) ' 
gh{l-h+3h-'-h^) ' 



PD- 1 
PD-2 
PD-3 
ED 

oo for the PD-1 case if u 7^ 0, 



Namely, the order in the correlation function survives as Nq 
whereas the order is destroyed in the other cases even for u^O. This result is just consistent with 
the numerical results for specific heat and the Polyakov lines studied in section 4. For A = case, 
the ordered- deconfinement phase is realized only in the PD-1 model at large g. 



5.2 The 2D XY model with temporal variables Uxo 

The discussion in the previous subsection is supported by considering another effective model that 
focuses on the dynamics of the temporal gauge field Uxo- To obtain it, we first replace the spatial 
gauge field Uxi by its average, 

Uxi^v (0 <«<!). (5.5) 

This replacement is supported for the PD-1 model by the result of the ID XY model (5.4) in the 
previous subsection. Then we fix the gauge to the temporal gauge by setting Uxo = 1 except for 
xo = No because of the periodic boundary condition. By keeping only the terms in the energy 
involving Uxj_,No,o, we have 

Ux±,No,0 = C^i-ifxJ, 

/TT ^ 
—^exp{A2nxY), 

X'_L ^ 

No 

A2DXY = gv'^^^^CrCOs{ipx^+i- iPxJ- (5.6) 
x± i=l,2r=l 

This is just the standard 2D XY spin model having the nearest-neighbor couplings among the XY 
spins cxp{iipx±)- Here the XY spin correlation function corresponds to the correlation function of 
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the Polyakov lines studied in section A. 



foj) = fp{x±)- 



(5.7) 



The spin stiffness in the effective 2D XY model is given as 



oo, 



PD 



1 



gv^Qi < 



2 2 
gv TT 

6 



PD 



2 



(5.8) 




PD 



ED 



3 



It is well known that the 2D XY model with finite spin stiffness has no long-range order, although 
the Kosterlitz-Thouless transition is possible. Only in the PD-1 case, the spin stiffness diverges for 
finite V and the model may have the order (Ux^^No.o) 7^ 0. This result (5.8) supports the procedure 
(5.1) for the ID XY model in the previous subsection for the PD-1 case. On the other hand, Eq.(5.4) 
supports the replacement (5.5) as we mentioned before. Thus studies of these two XY models give 
us the conclusion that only the PD-1 model at large g has the ordered (Coulomb) phase for the 
A = case. 

Studies of the two effective XY models in this section predict that the cases of PD-2, PD-3 and 
ED with A = have no Coulomb phase, as it was verified by the previous numerical calculations. 
Inclusion of A term, however, generates direct interactions between the spatial variables Uxi- For 
small A, Uxi's (i = 1, 2) can be integrated out perturbatively in powers of A. The obtained effective 
model contains nonlocal and multi-body interactions of (pxj_ 's which prefer the ferromagnetic order 
of (pxj_ ■ Therefore the effective model for A ^ 1 may not belong to the same universality class of the 
2DXY model with local interactions. Then the above result does not contradict the numerical result 
in the previous section, which shows the deconfinement phase exists in the PD-2 model with A = 1. 

6 Conclusion 

In this paper, we studied the nonlocal compact U(l) gauge theory on the 3D lattice, which "sim- 
ulates" gauge models coupled with massless/massive matter fields. The main contributions of the 
present paper may be the following two points: (i) MC simulations are feasible within reasonable 
computer time even for 3D lattice gauge theories with nonlocal interactions along one direction, 
and (ii) the measurements of E,C, Polyakov lines, Wilson loops, and local and nonlocal instantons 
give rise to clear and consistent results on the phase structure of the model. In particular, they 
distinguish the cases of ND(A = 1,0), PD-1(A = 1,0), PD-2(A = 1) with a CDPT from the other 
cases of PD-2(A = 0), PD-3,ED without CDPT in a definitive manner. 
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As explained in Sect. 2. 3, the results obtained in this paper are quite important for studies of the 
strongly-correlated electron systems like the high-Tc cuprates, the fractional quantum Hall effect, 
quantum spin models, etc. For example, in the t-J model of high-Tc superconductivity, by using 
the hopping expansion of holons and spinous at finite T with the continuous imaginary time, we 
derived an effective gauge theory, which is highly nonlocal in the temporal direction. The obtained 
effective theory has a similar action as Eq.(3.5) with Cr =constant and g oc n where n is the density 
of matter fields(holons and spinous) [4]. This corresponds to the ND model. Although the above 
effective gauge model is obtained for the system at finite T, we expect that a similar gauge model 
appears as an effective model at T = 0. The result that the ND model has a CDPT strongly suggests 
that the t-J model has the corresponding phase transition into the deconfinement phase, which is 
nothing but the charge-spin separated phase. 

In the deconfinement phase, all the three variables Ux^ (/U = 0,1,2) are stable, having small 
fiuctuations. The stability of the Lagrange multiplier Uxo means that the constraint (2.6) is not 
respected by the holons and spinous, the low-energy excitations of the system, whereas the stability 
of Uxi indicates that the gauge interaction between the holons and spinous can be treated pertur- 
batively. Therefore, quasi-particles in the CSS state are the holons, spinous and weakly interacting 
gauge bosons[4]. Of course, in order to present a definite "proof of the CSS, it is necessary to in- 
vestigate a gauge system with full isotropic nonlocal interaction, because the integration over holons 
and spinous generates nonlocal interactions not only in the temporal direction but also in the spatial 
directions and their combinations. 

Another interesting model related with the present one is the U(l) Higgs model coupled with 
the nonlocal gauge field. At present, it is believed that there is no phase transition in the 3D U(l) 
gauge-Higgs model with the ordinary local action if the Higgs field has the fundamental charge and 
its radial fiuctuations are suppressed [23]. However, the situation may be changed by nonlocal gauge 
interactions. The existence of the deconfinement phase in the present nonlocal gauge system without 
Higgs fields suggests that all the three phases, i.e., the confinement, Coulomb and Higgs phases, may 
be realized in the 3D nonlocal gauge model with a local coupling to a Higgs field. The deconfinement 
phase of the present model corresponds to the Coulomb phase. This problem is closely related with 
"doped holes" in the algebraic spin liquid which may be realized in certain antiferromagnetic spin 
models and materials in the spatial 2D lattice. We shall report on these problems in a separate 
publication [24] . 
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A High- and Low- Temperature Expansions 



In this appendix we study the behavior of the nonlocal gauge theory (3.5) in two regions of g by 
analytic methods; the region of small g{<< 1) by high-temperature expansion (HTE) in Sect. A. 1 
and the region of large g{>> 1) by low-temperature expansion (LTE) in Sect. A. 2. Once one obtains 
an approximate expression for the partition function Zq-, one can calculate various thermodynamic 
quantities like E and C in Eq.(4.1). 

A.l High-temperature expansion (HTE) for small g 

Let us consider the case of small g. Since the action Ar is proportional to g, one may expand the 
partition function Zr of (3.5) in powers of g as 

- {Arr 



Zr = J [dU] exp(Ar) = J [dU] ^ 
We obtain the expansion up to 0{g^) as 



n=0 



Zr^l+ {B2T + B2sX^) + B^Tg^ + {B^t + BiTS>? + BigX") g^ + 0{g''), 
where each coefficient is expressed as 

r "1 ^ 

But = ^[[dU]Y,''r{Vg + V,) ,q=ix,i,T), 

L q 

r -in 

LP 

^{Up + Up)-J2'^r{Vg + V,) 



BnS = 



B. 



4TS — 



4C2 

4! 



[dU] 



-\ 2 



LP q J 

To evaluate the above coefficients, we use the following formula for U(l) integral, 

J" dUxj^UxyJ^yv ~ ^xy^txv^mn- 

Then we obtain the following result; 

B2T = =2yQ2, Q2 = Y,cl, 

q r 

B2S = 5^1 = ^, 

p 

= X] X/ X/ X/ ^riCr2CT3Vx,i,Ti Vx+Tid,i,T2^':,i,r3^T3,Ti+r2 + C-C 

X i T\ T2 

N N N 



(A.l) 



(A.2) 



(A.3) 



(A.4) 



Tl = l T2 = l T3 = l 



28 



B4T = 



^4T 



fib 
-°4T 



'4T 



^ ^ ^ JJ^ ''T; • 4C2 ['J(jig2'^g3 94 + '^9193 '^13294 ~ '^91 92 '^92 93 '^93 94] 
' 91. 92, 93 ,94 t=l 



CriCr2Cr3Cr4(l ^rir3)^r4,ri+r2-r3Vi;,i,ri V^+T-i6,i,r2^'*)'^3K;+r36,i,T 



^ ^ 1 1 *^r3 '^T4 '^T4 , ti + T2 — ts ) 



4- 3- 2 

X,2 ri ,T2,T3 ,T4 

2F(Q4--02), 04- = 

Tl,r2,T3,r4 

y^ *^riCr2CT3Cr4Vi,i,riV^+ri6,i,T2^+(ri+T2)0,i,r3^'^)'^4^T4 + CC. 

ri .T2 .r3 .T4 

41^04+, ^4+ 



4! 
4! 



E 



*^Tl ^Tl *^T3 ^T4 ^T4 ,Tl +r2 +r3 5 



Tl,T2,T3,T4 



p 9 

-B4S = ^ ^ ^ 4C'2 [^pip2^P3P4 "I" ^P1P3^P2P4 ~ ^P1P2^P2P3'^P3P4] 



Pl,P2,P3,P4 

2 4 



(A.5) 



Some of them are illustrated in Fig. 14. 

Finally, we obtain the expression of the partition function and the free energy as follows, 

Zt = exp{-FV), 

F = -{2Q2 + X^)g^-4Qsg^+(^ + ^-AQi+-2QAg^ + 0{g'). (A.6) 
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In the following, we list up the values of Q2 and Q4 for A'' —> 00 





O2 




PD-1 


^ = 1.64493 


^ = 1.08232 


PD-2 


90 


9S0 = 1-00408 


PD-3 


= 1-01734 


1.00025 


ED 


0.15652 


0.018657 



(A.7) 



By using Zr of Eq.(A.6), we obtain the internal energy E and the specific heat C defined by Eq.(4.1) 
as 

E = -2{2Q2 + X^)g^-l2Q3g^ + A(^^ + ^-AQ4+-2Q4-^g^ + Oig'), 

C = 2(2Q2 + A^)^' + 24(535^ - 12 + ^ - 4Q4+ - 2Q4-) 5' + 0(5')- (A.8) 
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Figure 14: Illustration of the terms (A. 5) in HTE. Each rectangle represents Vx.i^r- The terms 
B^T, -84^ reflect the nonlocal nature; they are absent in models of local interaction characterized by 
= (r > 1). 




Figure 15: Comparison of HTE and the MC simulation. We compare the HTE results and the MC 
results of E and C for the PD-1 model with iV = 8 and A = 1. The curves indicated by En and Cn 
are the HTE results with the terms up to 0{g"). The dots arc the MC data. 
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Fig. 15 shows that, as one includes the higher-order terms, the HTE results approaches the MC 
result systematically. However, the approach is rather slow compared with the related models of 
local interactions like the 3D XY spin model or the 3D U(l) pure LGT. This is because the present 
nonlocal interactions generate various important higher-order terms in the HTE that are absent 
from the local models. 

Let us comment on the convergence of the HTE. As usual, the HTE is an expansion in the 
disordered (confinement) phase in which Ux^ fluctuates wildly. Equation (A. 6) shows that the 
convergence radius 5hte of the expansion is finite Qute 0) because both the harmonic numbers 
Q2 and Q4 appearing in the coefficients are finite. This means that there exists certainly the finite 
region < < 5hte confinement phase. We notice that if the long-range interaction is 

very strong such that Q2 = 00, the confinement phase may disappear. 



A. 2 Low- Temperature Expansion (LTE) for leirge g 



For large g, we evaluate Zr by the LTE. The LTE is an expansion in powers of g~^ around a fixed 
"lowest-energy configuration" of Uxfj, like J/^ju = 1, which gives rise to the global maximum of At- 
Let us expand as 

1 



= exp(i^^^) = 1 + iOx^ - -e' + o{etj, 



(A.9) 



where Ox^ is treated as 0{g ^/^) as we shall see. Aq- is expanded up to the second order in 9xfji in 
the following quadratic form; 

At = 4gQiV + 2gXV - gJ^Y. Ox^,Gx^,y.{\)ey, + O{0^), (A.IO) 

where the first two terms 4,gQ\V -\-2gXV come from the first term, unity, of R.H.S. of Eq.(A.9). Due 
to the gauge invariance, one may extend the region of 9xf_i from G (— 7r,7r) to 9^^ £ (—00,00) 
together with a gauge fixing. We take the temporal gauge 9xo = in the following calculation. Then 
we evaluate Zt by rescaling 9'^^ = g^^^9xi and performing Gaussian integration as 



Zt 



,(4Qi+2A)gy ■ 



^ r />oo /"OO "I / 

n / d^^i / d9 

^ V.J — 00 J — oo J V ^ ,, 



^ g(4Qi+2A)gV' ' 



r poc poo 1 / 

n / ^^xl / (^^'x2 exp -^^6';^G:rM,!/<^^yi. 

™ L J—oo J—oo J \ ^ n, 



(4Qi5 + 2Xg - lng)V - -TrlnG(A) 



= exp 

F = -(4Qi + 2A)5 + ln5 + 0(<?°). 



(A.11) 
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This gives 

E = _(4Qi + 2A)5 + 1 + 0(5-1), 

C = 1 + 0(5-'). (A.12) 

The higher-order terms in F are 0{g~'^) (n > 0) which may be calculated by the usual perturbation 
theory. 
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